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ABSTRACT 

This paper addresses the issue of how best to correct astronomical data for the wavelength-dependent effects of 
Galactic interstellar extinction. The main general features of extinction from the IR through the UV are reviewed, 
along with the nature of observed spatial variations. The enormous range of extinction properties found in the 
Galaxy, particularly in the UV spectral region, is illustrated. Fortunately, there are some tight constraints on the 
wavelength dependence of extinction and some general correlations between extinction curve shape and inter- 
stellar environment. These relationships provide some guidance for correcting data for the effects of extinction. 
Several strategies for dereddening are discussed along with estimates of the uncertainties inherent in each method. 
In the Appendix, a new derivation of the wavelength dependence of an average Galactic extinction curve from the 
IR through the UV is presented, along with a new estimate of how this extinction law varies with the parameter 
R = A(V )/E(B — V). These curves represent the true monochromatic wavelength dependence of extinction and, 
as such, are suitable for dereddening IR-UV spectrophotometric data of any resolution, and can be used to derive 
extinction relations for any photometry system. 

Subject headings: ISM:dust, extinction 



1. INTRODUCTION 

A precise knowledge of the wavelength dependence of inter- 
stellar extinction (i.e., the absorption and scattering of light by 
interstellar dust grains) and of any spatial variability in this de- 
pendence is important for two distinct reasons. First, the extinc- 
tion depends on the optical properties of the dust grains along 
a line-of-sight and potentially can reveal information about 
the composition and size distribution of the grains. Further, 
changes in the extinction from place to place may reveal the de- 
gree and nature of dust grain processing occurring in the ISM. 
Second, the wavelength dependence of extinction is required 
to remove the effects of dust obscuration from observed en- 
ergy distributions, since most astronomical objects are viewed 
through at least some small amount of interstellar dust. Spa- 
tial variations in the extinction potentially limit the accuracy to 
which energy distributions can be "dereddened." Such uncer- 
tainties might be acceptably small for very lightly reddened ob- 
jects, but quickly can become debilitating along modestly red- 
dened sightlines. 

This paper follows from a talk given at the meeting "Ultravi- 
olet Astrophysics Beyond the IUE Final Archive" (Fitzpatrick 
1998) and addresses the second point raised above, i.e., the 
correction of energy distributions for the effects of interstel- 
lar extinction. The goals are to provide a summary of what is 
known about the wavelength dependence and spatial variability 
of interstellar extinction and to present strategies for the best re- 
moval of the effects of extinction from astronomical data. In §2, 
I briefly review the chief features of IR-through-UV extinction 
and the nature of the known spatial variations. These variations 
are illustrated with a set of optical/UV extinction curves which 
span the (currently) known extreme limits of extinction varia- 
tions. Constraints on the wavelength dependence of extinction 
and some general correlations between extinction curve shape 
and interstellar environment are also noted in §2. Strategies for 
dereddening are discussed in §3 along with estimates of the ap- 
propriate uncertainties which ought to be incorporated into an 
error analysis. Some final comments regarding extinction along 



extragalactic sightlines and at wavelengths shortward of 1 150 A 
are given in §4. 

The Appendix of this paper describes the construction of 
a new estimate of the shape of a mean and an environment- 
dependent IR-through-UV extinction curve. These curves, 
which may be obtained from the author, are suitable for dered- 
dening multiwavelength spectrophotometric data, such as has 
become available with the Hubble Telescope's Faint Object 
Spectrograph, and for deriving extinction relationships for any 
photometric system. 

2. SPATIAL VARIATIONS IN GALACTIC EXTINCTION 
2.1. Variations Abound ... 

Figure 1 shows several estimates of the shape of a "mean" 
Galactic extinction curve from the far-IR through the UV re- 
gion. The curves are presented in a commonly used normal- 
ization scheme, E(X-V)/E(B-V), and plotted against inverse 
wavelength. The significance of the different curves will be dis- 
cussed later in this paper. For the present, these curves serve to 
illustrate the overall characteristics of interstellar dust extinc- 
tion. Typically, the extinction rises through the IR with a power 
law-like dependence (see Appendix), rolls over slightly in the 
optical region ("knee"), shows a prominent feature at 2175 A 
in the near-UV ("bump"), and has a sometimes steep rise in the 
far-UV ("fuv rise"). Figure 1 does not attempt to show the rela- 
tively narrow optical/IR features which may be associated with 
interstellar dust, such as the diffuse interstellar bands in the op- 
tical or the solid state absorption (or emission) features in the 
IR (e.g., the 9.8 /im silicate feature). 

IR and optical extinction traditionally have been studied us- 
ing ground based broad- or intermediate- band photometry and 
sightline-to-sightline variations have been recognized for a long 
time (e.g., Johnson 1965). These are often characterized using 
the parameter/? = A(V)/E(B-V), i.e., the ratio of total to se- 
lective extinction at V. The value of R ranges between about 2.2 
and 5.8 for sightlines along which UV extinction has also been 
measured, with a mean value of ~3.1 for the diffuse interstellar 
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medium. It is believed that the shape of the IR extinction law 
longward of ^8000 A (the "power law" region) may be invari- 
ant and that the observed range in R is due to spatial variations 
in the steepness of the extinction in the optical region around 
the "knee" (e.g., Cardelli, Clayton, & Mathis 1989; Martin & 
Whittet 1990). 

Measurements of UV interstellar extinction have relied on 
satellite or rocket spectrophotometric data, and date back to 
the original discovery of the 2175 A bump nearly 30 years ago 
(Stecher 1969; Bless & Savage 1970). These studies immedi- 
ately revealed that sightline-to-sightline differences exist within 
the Galaxy in the detailed properties of UV extinction curves. 
Important results on the degree of spatial variability of UV ex- 
tinction were subsequently derived from the large photomet- 
ric databases accumulated by the TD-1 and ANS satellites (e.g. 
Koornneef 1978; Kester 1981; Meyer & Savage 1981). 

The most graphic illustrations of spatial variations in UV ex- 
tinction were provided by the IUE satellite. With its complete 
spectral coverage between 1 150 and 3200 A and relatively high 
spectral resolution (^6 A in the low-resolution mode) IUE was 
well-suited for the study of UV extinction and many studies 
have utilized IUE data to determine the wavelength dependence 
of extinction along one or more Milky Way lines-of-sight (e.g. 
Bohlin & Savage 1980; Witt, Bohlin, & Stecher 1980; Seab, 
Snow, & Joseph 1981; Morgan, McLachlan, & Nandy 1982; 
Hecht et al. 1982; Massa, Savage, & Fitzpatrick 1983; Witt, 
Bohlin, & Stecher 1984; Massa & Savage 1984; Franco, Mag- 
azzii, & Stalio 1985; Aiello et al. 1988; and many others). Fig- 
ure 2 shows analytical representations for ~80 UV extinction 
curves derived from IUE data (solid curves), illustrating the 
wide range of properties observed for Galactic UV extinction. 
These curves are taken from the catalog of Fitzpatrick & Massa 
(1990), with a few additions, including the lines-of-sight toward 
HD 210121 (Welty & Fowler 1992) and HD 62542 (Cardelli & 
Savage 1988). The basis for the analytical representation will 
be discussed further below. 

The overall impression that has come from the last 30 years 
of extinction studies is that there is a bewildering variety of 
IR-through-UV extinction curves and that what has often been 
called "peculiar" extinction is actually quite common. This 
bodes poorly for the prospects of precisely dereddening astro- 
nomical energy distributions. In particular, the range of extinc- 
tion properties seen in the UV (Figure 2), where the effects 
of extinction are strongest, may represent a bonanza of infor- 
mation for those who study the properties of interstellar dust 
grains, but is a potential disaster for those whose main interest 
in UV extinction is the removal of its effects. 

2.2. ... But There is Some Order Amid the Chaos 

The recognition that interstellar extinction is spatially highly 
variable is critically important for both the study of dust grains 
and for the purposes of dereddening. Fortunately for the latter, 
it has been realized that within the seemingly endless variety of 
observed extinction curves there are actually some constraints 
on the wavelength dependence and some links between extinc- 
tion in the UV, optical, and IR. 

This began with Savage (1975) who showed, using OAO-2 
satellite data, that the shape of 2175 A bump can be reproduced 
well with a Lorentzian profile. This result — which provides 
important clues about the nature of the 2175 A feature — al- 
lows a functional representation of the bump, and permitted 
Seaton (1979) to present an estimate of the average Galactic 



UV extinction curve as a simple analytical formula. This curve 
is shown in Figure 1 by the dashed curve, extending from 2.7 
to 8.7 /an -1 . 

Following the lead of Savage and Seaton, Fitzpatrick & 
Massa (1986, 1988, 1990; collectively referred to hereafter as 
FM) examined IUE extinction curves for many lines-of-sight 
and found that all these curves could be fitted extremely well 
by a single analytical expression with a small number of free 
parameters. This expression consists of (1) a Lorentzian-like 
bump term (requiring three parameters, corresponding to bump 
width 7, position xO, and strength c3), (2) a far-UV curvature 
term (one parameter c4), and (3) a linear term underlying the 
bump and the far-UV (two parameters cl and c2). This set of 
basis functions is shown in Figure 3. With the proper choice 
of the 6 parameters, essentially all UV extinction curves can be 
reproduced to within the uncertainties inherent in the data. For 
most sightlines, this number of free parameters can be reduced 
to 4, since it has been seen that the position of the 2175 A bump 
is nearly invariant and that the intercept and slope of the linear 
term are tightly correlated (FM; see also Carnochan 1986 and 
Jenniskens & Greenberg 1993), such that the linear term ap- 
pears to pivot about the point [1/A,f^j] ~ [3,2]. The UV 

portions of the curves shown in Figure 2 (A < 2700 A) were all 
constructed using these functions. There is, thus far, no con- 
vincing evidence for substructure in UV extinction curves. 

The FM results show that some order underlies the wide 
spectrum of wavelength dependences seen in the UV, but do 
not necessarily permit a more accurate dereddening of an ob- 
ject along any given line of sight. Fortunately it has been ob- 
served that extinction curves do not vary randomly across the 
sky, but show distinct regional signatures (e.g., Meyer & Sav- 
age 1981; Morgan, McLachlan, & Nandy 1982; Panek 1983; 
Massa & Savage 1984; Clayton & Fitzpatrick 1987) and it is 
clear that extinction properties reflect the physical conditions 
and past processing histories of the environments in which the 
dust grains reside. 

The work of Cardelli, Clayton, & Mathis (1988, 1989; here- 
after CCM) provides a link between one measure of dust grain 
environment and the wavelength dependence of U V-through-IR 
extinction. CCM found that the shape of UV extinction curves 
correlates with the parameter R (defined in §2.1 above). This 
suggests that — although there is considerable scatter — ex- 
tinction curves from the UV through the IR can be character- 
ized as approximately a one-parameter family dependent on R. 
Thus, if the value of R can be determined (from optical and IR 
photometry), then the properties of the the entire UV-IR ex- 
tinction curve can be predicted. The essence of the CCM result 
is illustrated in Figure 4. Four representative extinction curves 
are shown, each determined by the value of R listed at the right- 
hand side of the curves. This figure shows that extinction curves 
which are seen to be very "flat," or "grey," in the UV roll over 
strongly in the optical and are characterized by large values of 
R. Steep UV curves remain steep in the optical and are char- 
acterized by small values of R. The CCM curve for R = 3.1 is 
shown in Figure 1 by the dotted line. 

The results of CCM are important in 3 ways: (1) They indi- 
cate that some of the spatial variations seen in extinction curves 
behave coherently and systematically over a wide wavelength 
range, thus potentially allowing for the consistent dereddening 
of multiwavelength energy distributions. (2) The observed de- 
pendence on R allows for the definition of a meaningful aver- 
age Galactic extinction curve. Existing datasets of extinction 
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curves, particularly for the UV, are biased toward "interesting" 
regions where extinction properties are extreme. Computing a 
simple mean curve from these data does not necessarily yield 
a reasonable estimate of a Galactic mean. Since it is well- 
established, however, that the appropriate mean value of R for 
the diffuse ISM is -3.1 (e.g., Schultz & Wiemer 1975; Whit- 
tet & van Breda 1980; Rieke & Lebofsky 1985), a reasonable 
definition of a mean Galactic extinction law is that which corre- 
sponds to the case where R = 3.1. (3) They demonstrate a gen- 
eral correlation between dust grain environment and the wave- 
length dependence of extinction, since large values of R are 
generally found in dense environments where dust grain growth 
is thought to occur. The "greyness" of the UV/optical extinc- 
tion curves for large R is consistent with a larger than normal 
dust grain population. These results have an important bearing 
on how to correct for the effects of extinction, discussed in the 
following section. 

3. CORRECTING FOR THE EFFECTS OF EXTINCTION 

Given the observed complexity of Galactic UV extinction, 
what is the best way to deredden an observed energy distribu- 
tion? There are three different possibilities, which involve the 
use of a global mean extinction curve, an /^-dependent curve, 
and a sightline-specific curve. These three alternatives and their 
inherent uncertainties are discussed below, and a new derivation 
of the shape of the average IR-through-UV extinction curve is 
presented. A more thorough discussion of the uncertainties in- 
volved in UV extinction corrections is given by Massa (1987). 

3.1. Using a Global Mean Curve 

If there is no specific information available about the wave- 
length dependence of the extinction along a sightline of interest 
then the only alternative is to adopt some globally defined mean 
curve and perform a realistic error analysis. This is the least at- 
tractive of the three cases discussed here, but is by far the most 
commonly used. The average Galactic extinction curves from 
either Seaton (1979; see Figure 1) or Savage & Mathis (1979) 
are often adopted for dereddening UV data. These should now 
be superseded by the use of an /^-dependent curve computed 
for the case R = 3.1 (see §2.2). Figure 1 shows the R = 3.1 
curve from CCM (dotted line) and a new determination of the 
R = 3.1 case (solid line) which is described in the Appendix to 
this paper. This new curve aims to reproduce the detailed wave- 
length dependence of the R = 3.1 extinction law and has been 
constructed to account for bandpass effects properly in opti- 
cal/IR extinction data and to reproduce the observed broad-, 
intermediate-, and narrow- band extinction measurements. The 
curve is thus suitable for dereddening multiwavelength spec- 
trophotometric observations and can be used to derive the aver- 
age extinction relationships in any photometric system. 

A complete evaluation of the likely error in a dereddened rel- 
ative energy distribution m(A - V) requires the propagation of 
the uncertainty in E(B-V) and the often-neglected but often- 
dominant uncertainty in the adopted mean curve. This is given 
by 

<a-v) = 4 ( b-v) x k(X-V) 2 +E(B-V) 2 x aj {X _ v) (1) 

where k(X-V) represents the normalized extinction curve E(\— 
V)/E(B-V). If absolute fluxes are desired, then the uncertainty 
in the assumed value of R must also be incorporated (two addi- 
tional terms on the righthand side of eq. 1, similar to the current 
terms but with R substituted for k(X-V)). 



The large database of ANS satellite extinction measurements 
from Savage et al. (1985) shows that the 1-er scatter at 1500 A is 
<7> (1 5_ y) = 0.74, based on -400 sightlines with E(B-V) > 0.5. 
(Only relatively large values of E(B-V) were considered in 
order to minimize the effects of spectral mismatch error and 
random noise on the a measurement.) To extend this analy- 
sis to other wavelengths, we computed the standard deviation 
at each wavelength point for the 80 curves shown in Figure 
2, and then scaled the result to match the ANS value at 1500 
A. (The actual standard deviation at 1500 A for the 80 IUE 
curves is somewhat higher than the ANS result because of the 
bias in Figure 2 toward extreme extinction curves). The resul- 
tant values of <Jk{\-v) are shown by the thick dotted curve near 
the bottom of Figure 2 (labeled "a") and are listed at selected 
wavelengths in the third column of Table 1 . This estimate of 
crk(\-v) should be adopted whenever the average Galactic ex- 
tinction curve is used for dereddening an observed energy dis- 
tribution. The uncertainties approach zero for 1/A < 3 /im" 1 
due to the curve normalization. The quantity E(B-V) is usu- 
ally derived directly from photometry and often has an easily 
quantifiable uncertainty; thus the computation of <Tm(\-v) fr° m 
eq. 1 is straightforward. 

Sometimes E(B - V) is not known a priori and is estimated 
by "ironing out" the 2175 A extinction bump using an assumed 
extinction curve shape (e.g., see Massa, Savage, & Fitzpatrick 
1983). In many cases, uncertainties in E(B-V) as small as 0.01- 
0.02 mag for moderately reddened objects have been quoted 
from this process. This is incorrect! The normalized "height" of 
the 2175 A extinction bump — which is the quantity that is im- 
portant in the ironing-out process — has a \-a scatter of about 
±20% around its mean value (from the data in Figure 2 and 
Savage et al. 1985). Thus the uncertainty in an E(B-V) mea- 
surement derived from the bump must be considered to have 
a similar relative uncertainty. Since bump strength does not 
correlate well with other aspects of UV extinction, such as the 
slope of the linear component or the strength of the far-UV rise 
(FM), the uncertainty in an energy distribution dereddened this 
way can be estimated by using eq. 1 with (Je(b-v) — 20% and 
the values of (Jk(\-v) from Table 1 . The result of this calcula- 
tion needs to be modified somewhat by removing the signature 
of the 2175 A bump. The final estimate of cr^ l(A _ l/) computed 
this way is listed in the last column of Table 1 . Multiplying 
these values by E(B - V) gives the uncertainty in a dereddened 
relative energy distribution. The values in Table 1 agree well 
with estimates a 2 m{X _ v) at 1500, 1800, 2200, 2500, and 3300 

A derived from the ANS data of Savage et al (1985), and are 
consistent with the analysis of Massa (1987). 



3.2. Using a R-Dependent Mean Curve 

If the value of R has been measured along a sightline of in- 
terest, or if some general information is available about the dust 
grain environment (e.g., dense region or diffuse region?), then 
the uncertainties in dereddening can be reduced somewhat by 
adopting an /^-dependent extinction curve computed at the ap- 
propriate value of R. For the original 29 sightlines used by 
CCM, the standard deviation of the individual observed extinc- 
tion curves minus the computed /^-dependent curves is about a 
factor of 2 smaller than when the individual curves are com- 
pared to the average Galactic (R = 3.1) curve. Thus the un- 
certainties in a dereddened energy distribution can be estimated 
using eq. 1 and adopting 0.5 times the (Jk(\-v) values from Table 
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1 . If R is poorly measured or only estimated from environmen- 
tal factors, then allowance must be made for this in the error 
analysis. 

Using an appropriate /^-dependent curve can slightly improve 
the accuracy of results from ironing out the 2175 A bump. The 
values of cr^ (A _ v) can be estimated as above in §3.1, using the 
same value of (Je(b-v) — 20% and 0.5 x a^x-v) from Table 1 . 
The resultant values of <y^x-v) are arj out 85% as large as those 
listed in the last column of Table 1 . The gain from using the 
/^-dependent curves is only marginal because the normalized 
bump heights are highly variable and not correlated with R. 

It is important to note that, even when R and E(B - V) are 
well-determined, the use of an /^-dependent extinction curve 
only reduces — but does not eliminate — the wavelength de- 
pendent dereddening error. The standard deviation of the ob- 
served vs. /^-dependent curves does not go to zero, even for 
the sample of 29 sightlines used by CCM to define the R- 
dependence, because the UV/optical curves are not really a one- 
parameter family dependent only on R. Curves derived for dif- 
ferent sightlines with the same values of R show a wide range 
of properties, including differences in the strengths of the bump 
and the far-UV rise. Thus, as noted by CCM, the CCM formula 
reproduces a general trend, but does not provide particularly 
good fits to individual extinction curves, even when the value 
of R is well-determined. 

Apart from this intrinsic scatter around the mean R relation, 
the accuracy of the CCM results is limited in the IR/optical re- 
gion due to bandpass effects with the broadband Johnson fil- 
ters used to measure the extinction. As a result, CCM tends 
to overestimate the level of extinction in the near-IR and blue- 
visible. In the Appendix, a new derivation of the /^-dependence 
of IR-through-UV extinction is presented. This corrects for the 
systematic bandpass effects in CCM, but is still plagued by the 
same uncertainties caused by the "cosmic" scatter of extinction 
properties around the mean R relation. 

3.3. Using a Sightline-Specific Curve 

Without question, the best of all possible dereddening 
choices is to use an extinction curve derived explicitly for the 
sightline of interest. This advice is not so unhelpful as it may 
sound, because the archives of the IUE satellite may well con- 
tain normal stars in the direction (and at the distance) of the 
object of interest, from which extinction curves may be derived 
(e.g., halo stars in the case of extragalactic sightlines). The ac- 
curacy of the dereddened energy distribution may thus be lim- 
ited chiefly by the uncertainties in the "pair method" of con- 
structing the extinction curve (see, e.g., Massa et al. 1983) and 
by the determination of E(B-V) for the object of interest. For a 
well-determined curve, "ironing out" the bump — in the cases 
when E(B-V) is not known explicitly — should yield the same 
degree of accuracy as for the cases when E(B-V) is known a 
priori. An example of such an approach is the determination 



of the intrinsic energy distributions of rare or exotic objects lo- 
cated in clusters with more mundane objects. Massa & Savage 
(1985) utilize this technique, using extinction curves derived 
from main sequence B stars in open clusters to determine the 
continuum shapes of cluster O stars. 

4. FINAL COMMENTS 

The results discussed above strictly pertain to Milky Way ex- 
tinction for wavelengths long ward of about 1150 A (the IUE 
satellite cutoff). Observers of extragalactic objects must con- 
tend with foreground Milky Way extinction as well as extinc- 
tion due to dust grains local to the object of interest. It is rea- 
sonable to assume that the properties of interstellar extinction 
in external galaxies are fully as complex as found for the Milky 
Way. However, only for the nearby Large Magellanic Cloud 
(LMC) do observed extinction curves for individual sightlines 
approach the accuracy needed to study spatial extinction vari- 
ations (see Fitzpatrick 1998 and references within). The prop- 
erties of Milky Way extinction along the high galactic latitude 
sightlines relevant to most extragalactic observations are also 
not well-determined. One of the few quantitative studies of 
high latitude extinction, that by Kiszkurno-Koziej & Lequeux 
(1987) utilizing ANS satellite data, suggests that halo extinction 
is slightly steeper in the far-UV and has a weaker 2175 A bump 
than the typical extinction found in the plane of the Galaxy. At- 
tempts to eliminate galactic foreground reddening by ironing 
out the bump may lead to systematic underestimates of the total 
amount of foreground extinction. 

The properties of extinction in the wavelength range 912- 
1150 A are poorly known, but the work that has been done 
(e.g., York et al. 1973; Snow, Allen, & Polidan 1990) suggests 
that the far-UV rise continues at least down to wavelengths of 
950 A or so. An intrinsic problem with measuring extinction in 
this wavelength range is the high density of H I and H2 absorp- 
tion lines which mimic continuous absorption in low-resolution 
data. CCM suggest that the wavelength dependence of extinc- 
tion in this region is reasonably well-represented by an extrapo- 
lation of the fitting functions used in the IUE wavelength range. 
At the present, the best dereddening strategy for 912-1150 A 
data would be to adopt such an extrapolation, based if possible 
on IUE extinction curves derived for the sightline of interest. 
The uncertainties listed in Table 1 for 1000 A are based on this 
simple extrapolation. 

I thank my partner-in-interstellar-grime, Derek Massa, for 
many years of enjoyable collaboration. I also acknowledge the 
contributions of the late Jason Cardelli to the study of inter- 
stellar extinction. Future progress in the field will proceed at a 
slower pace without Jason's enthusiasm and scientific insight. 
This work was supported in part by NASA grant no. NAG5- 
7113. 



APPENDIX 

DERIVING A /^-DEPENDENT IR-THROUGH-UV EXTINCTION CURVE 

In this Appendix a new estimate of the wavelength dependence of Milky Way extinction from the IR through the UV is derived. The 
motivation for this is the need for a multiwavelength extinction curve suitable for dereddening IR-through-UV spectrophotometric 
data now available from, for example, the Hubble Space Telescope's Faint Object Spectrograph (see, e.g., Fitzpatrick & Massa 1998 
and Guinan et al. 1998). The specific goals are 1) to produce a curve for the R = 3.1 case which reproduces existing narrow band 
extinction measurements and the principal Johnson and Stromgren photometric extinction relations and 2) to define the /^-dependence 
of the extinction law, relying principally on spectral regions where this dependence is well-determined. Existing estimates of the 
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average wavelength dependence of optical/near-IR extinction do not accurately portray the true shape of the curve due to bandpass 
effects arising from the broadband photometry used to derive the curves. This deficiency will be corrected by using synthetic 
photometry. 

The UV Region 

CCM noted a correlation between the value of R~ l (as determined from IR/optical photometry) and the values of A(X)/A(V) (i.e., 
total extinction at A normalized by total extinction at V) at UV wavelengths using parametrized extinction curves for a subset of 
29 stars from the FM sample. They consequently derived a complex polynomial expression to reproduce the wavelength- and R- 
dependences. We approach the problem slightly differently and note that the essence of the CCM result is two correlations involving 
the coefficients of the FM fitting function originally used to parametrize the extinction curves studied by CCM. We thus derive the 
functional forms of these correlations and retain the FM fitting function to compute the values of E(X-V)/E(B-V) for wavelengths 
A < 2700 A. 

Figure 5 illustrates the two correlations. The top panel shows a plot of the FM parameter c2 (representing the slope of the linear UV 
extinction component) vs. R~ l (filled circles) along with a least squares estimate of the linear relationship between the two (dashed 
line). The data are for a subset of the FM sample for which IR photometry could be used to deduce the value of R (see below), 
with the addition of the sightline toward HD 210121 (Welty & Fowler 1992; Larson, Whittet, & Hough 1996). This represents the 
only convincing correlation between R and the properties of the UV extinction curve (see, e.g., Jenniskens & Greenberg 1993). The 
equation of the least-squares fit to the data is 

c2 = -0.824 + 4.717 x/T 1 (Al) 

No physical interpretation or significance is placed on the form of the adopted functional relationship between R and c2; it is merely 
that which best reproduces the observed correlation. We note, however, that the intrinsic slope of extinction in the optical region is 
proportional to R~\ and that eq. Al thus simply implies that the slopes of the extinction curves in the UV and optical vary together 
in a linear manner — as optical extinction steepens, UV extinction steepens. 

The bottom panel of Figure 5 shows the well-known relationship between slopes (c2) and intercepts (cl) of the linear background 
component (FM; Carnochan 1986; Jenniskens & Greenberg 1993). The data are for the full set of 80 curves from FM, plus HD 
210121, and minus the Orion Nebula stars (HD 36982, 37022, 37023, and 37061) which suffer from scattered light contamination 
(which mainly affects the linear intercept). The equation of the least squares fit to this relationship (dashed line), is 

cl =2.030-3.007 xc2 (A2) 

We use the following mean values for the other four parameters required to specify UV extinction curves with the FM formula: xQ 
(bump position) = 4.596 /xm" 1 ; 7 (bump width) = 0.99 ^m" 1 ; c3 (bump strength) = 3.23; and c4 (FUV curvature) = 0.41. 

The Optical/IR Region 

A common way of representing the wavelength dependence of extinction in the optical and IR regions is to draw a smooth curve 
through the values of E(X-V)/E(B-V) (or some other curve normalization) derived from Johnson UBVRIJHKLM photometry, 
assuming that these represent the monochromatic values of the extinction at the filter effective wavelengths A e // (e.g., CCM, Martin 
& Whittet 1990). This procedure does not accurately yield the monochromatic extinction curve because it ignores the wavelength 
dependence of extinction across the width of an individual filter. Extinction decreases towards longer wavelengths in the optical/IR 
region and therefore the value of A e // for a reddened star is shifted to a longer wavelength than for an identical unreddened star. An 
extinction measurement made by comparing photometric indices for two such stars (i.e., the "pair method") will always overestimate 
the monochromatic extinction in the neighborhood of A e // because stellar emergent fluxes in the optical/IR decrease towards longer 
wavelengths (for early-type stars). The magnitude of this effect depends on how much the true extinction curve varies across the 
filter. A broadband E(X-V)/E(B-V) measurement will also be influenced by the intrinsic energy distributions of the stars used in 
the pair method and by the total amount of extinction. 

The approach taken here is to find the wavelength dependence of the optical/IR extinction curve which reproduces the photometric 
extinction measurements when synthetic photometry of an artificially reddened stellar energy distribution is compared with that of 
identical but unreddened star. To represent the stellar energy distribution, an ATLAS9 model atmosphere from R.L. Kurucz (with 
T e ff = 30000 K and logg = 4.0) is used, along with an adopted value of E(B-V) = 0.5. Synthetic photometry is performed to yield 
measurements in the broadband Johnson UBVRIJHKLM system and the intermediate band Stromgren uvby system. The first step is 
to determine the curve shape for the mean R = 3.1 case and then to define the nature of the variation with R. The first two columns 
of Table 2 list the photometric extinction ratios (in both the Johnson and Stromgren systems) and their observed values, which the 
R = 3.1 extinction curve is required to reproduce. References to the observations are given in Column 3. In addition, the curve is 
constrained to reproduce the mean narrowband measurements in the 3400-7900 A region published by Bastiaansen (1992), which 
are assumed to represent the monochromatic values of the extinction for the case R = 3. 1 at the filter central wavelengths. 

Figure 6 shows the shape of the monochromatic extinction curve which best satisfies these requirements (thick solid curve), plotted 
as total extinction A(A) normalized by E(B-V). The arbitrarily scaled profiles of the Johnson and Stromgren filters are indicated, 
and the Bastiaansen data shown by the small plus signs ("+"). The agreement between the R = 3. 1 curve and the Bastiaansen data is 
clear from the figure. The fourth column of Table 4 gives the values of the photometric extinction ratios produced by the curve, also 
in very good agreement with the observations. 

For wavelengths longward of 2700 A(l/A < 3.7pm~ 1 ), the R = 3.1 curve is defined by a cubic spline interpolation between a set 
of optical/IR anchor points (filled circles) and a pair of UV anchor points (filled squares). The values of the UV anchors (at 2700 A 
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and 2600 A) are determined by the FM fitting function for the case R = 3.1 (see above) and assure a smooth junction between the 
optical and UV regions at 2700 A. The anchor at /im" 1 is fixed at (i.e., no extinction at infinite wavelength) and the values of the 
other six optical/IR points were adjusted iteratively to find the curve shape which best reproduced the extinction observations. The 
wavelengths chosen for the spline anchor points are somewhat arbitrary, although points in the IR, the optical normalization region, 
and the near-UV are clearly required. The wavelengths and values of the spline anchors for the R = 3.1 curve are given in Table 3. 

Note the slope of the derived extinction curve does not approach zero as 1 /A approaches zero. We sacrificed this physically 
expected requirement in order to achieve a better fit to the Johnson IR photometry and to preserve the simplicity of the fitting 
procedure. A zero slope could have been guaranteed by, for example, adopting a power law to represent the extinction; but no 
single power law can reproduce the IR photometry to an acceptable level. The new curve should be treated as very approximate at 
wavelengths beyond the limit of the M band (i.e., at A > 6/Ltm). 

The overall /^-dependence of the optical/IR curve is relatively easy to incorporate by the following adjustments in the spline 
anchors: (1) the UV points are computed by the FM fitting function using the coefficient values given above, including the R- 
dependent value of c2; (2) the IR points at 1/A < 1/im" 1 are simply scaled by R/3.1, since the shape of the far-IR extinction is 
believed to be invariant (see §2.1 and below); and (3) the optical points are vertically offset by an amount R — 3.1, with slight 
corrections made to preserve the normalization. Without these corrections (which are less than 0.015 over the range R — 2 to 6) 
the extinction curves would drift away from the standard normalization, i.e., E(B-V) = 1, by a few hundredths of a magnitude as 
R departed from the value 3.1. Table 4 gives formulae for computing the /^-dependent values of the optical spline anchors. The 
corrections just noted are manifested in the departure of the linear term from a value of 1.0 and in the higher order term for the 41 10 
A point. 

It is not obvious a priori that the above adjustments in the spline should produce the proper curve shapes in the relatively large 
gaps between the 4 points in the optical normalization region and the UV or IR regions. However, it will be shown below that the 
resultant curves are in good agreement with observations. 

The Full R-Dependent Curve 

Figure 7 shows the full wavelength range of the IR-through-U V /^-dependent extinction curves derived here, for four representative 
values of R (thick solid and dashed lines). For A < 2700 A the curves are computed using the FM fitting function and for A > 2700 
A the curves are spline interpolations between the /^-dependent spline anchors. The CCM results are shown for comparison by the 
thin dotted lines for the same four R values. In the UV region, the new results and the CCM curves are similar for R < ^4.0, but 
diverge for larger/?. The main reason for this discrepancy lies in the UV linear extinction component. While CCM adopted a linear 
relationship between the slope c2 and intercept cl, their relation — shown by the dotted line in the bottom panel of Figure 5 — does 
not agree well with the observations. The large discrepancy at small values of c2 (i.e., large R) produces the difference seen in Figure 
7. 

A closeup comparison between the new results and CCM in the optical/IR region for the R = 3 . 1 case can be seen in Figure 6, where 
the CCM curve is indicated by the dotted line. The disagreement in the region near the Johnson R filter illustrates the bandwidth 
effect discussed above. The CCM curve is fixed at A(X)/E(B-V) = 2.32 at 7000 A, which is taken as A c // for the Johnson R. The 
new curve is constructed to yield the same value of A( A) /E(B- V) but for synthetic photometry with the Johnson R. The discrepancy 
is particularly large for this filter because of the steep change in extinction across its relatively broad profile. This effect also accounts 
for the lesser discrepancies in the regions of the U, B, I, and J filters. O'Donnell (1994) presented a revision to the CCM formula 
in the region of the U, B, and V filters based on the Stromgren photometric indices. This revision uses the same stategy as CCM of 
fixing the extinction ratios at the filter effective wavelengths, but more closely resembles the newly derived monochromatic extinction 
curve because the bandwidth effects in the narrower Stromgren filters are smaller. 

Note that the technique used here to produce /^-dependent curves can be used to construct "customized" UV/optical extinc- 
tion curves for sightlines with well-defined UV extinction properties. It is simply necessary to substitute the measured values of 
A(X)/E(B-V) at 2700 A and 2600 A for the UV spline anchors listed in Table 3 and then perform the cubic spline interpolation to 
determine the optical portion of the curve, which smoothly joins the measured UV curve at 2700 A. This method was used in Figure 
2 to extend the IUE curves into the optical region. 

The assertion that the shape of the IR extinction law is invariant can be tested by comparing results obtained from the /^-dependent 
curves in Figure 7 — which were constructed based on this assumption — with observations. This comparison is made in Table 
5 for a number of sightlines which span the observed range in R. The "model" data listed for each sightline are the results from 
synthetic photometry on model stellar energy distributions artificially reddened using the /^-dependent curves, with E(B-V) values 
from column 2 of the table and R values from column 7. The listed values of R are those that best reproduce the observed color 
excesses and the uncertainties result from assuming uncertainties of ±0.02 mag in E(B-V) and ±0.05 mag in the IR E(V- A). The 
agreement between observations and model values is excellent and there are no strong systematic trends evident. The results are thus 
consistent with an invariant IR extinction curve. At wavelengths greater than ~1 /im, the extinction curve roughly resembles a power 
law with an index of ~1.5. This is similar to that adopted by CCM (1.6), but much flatter than that of Martin & Whittet (1990; ~1.8). 

Values of R are sometimes estimated from the relation R ~ 1.1 x E(V-K)/E(B-V). Exact relationships between R and the IR 
color excesses can be derived for the new /^-dependent curves. These are given by 

E(V-J) 

R = 1 .39 x — - 0.02 (A3) 

E(B-V) 

R= 1.19 x ^— — — Uo.04 (A4) 
E(B-V) 
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R=l.l2x 



#=1.07 x 



E(V-K) 
E(B-V) 
E(V-L) 
E(B-V) 



-0.01 



+ 0.02 



(A5) 



(A6) 



The coefficients in each of these equations actually depend on the value of E(B-V) itself but, over the range E(B — V) = - 2.0, vary 
by only several hundredths. The values in the equations are the results for E(B-V) = 0.5. 

At the blue end of the optical region, the interpolation between the optical and UV spline anchor points can be tested by compar- 
ing predictions of the extinction indices E(U-B)/E(B-V) and E(cl)/E(b-y) with photometric measurements. Figure 8 shows this 
comparison, with the difference between the observed values and predicted values — derived from synthetic photometry of "cus- 
tomized" extinction curves from the FM sample — plotted against the slope of the linear component c2. The curves were produced 
as described above, by adopting the observed values of A(X)/E(B-V) at 2600 A and 2700 A as the UV anchor points. The random 
scatter (measurement noise) in both indices is large, but no systematic trends are seen and the model curves appear to reproduce 
at least the integrated properties of extinction well in the near-UV region. Over the range of c2 values shown, the model values of 
E(U-B)/E(B-V) range from about 0.6 to 0.8 and the values of E(cl)/E(b-y) from about 0.0 to 0.3. 

In summary, the /^-dependent curves derived here reproduce existing photometric and spectrophotometric measurements and 
provide a good estimate of the true monochromatic wavelength dependence of interstellar extinction in the IR-through-UV regions. 
These curves should be preferred for dereddening UV, optical, and near-IR spectrophotometry. Since the new curves give the detailed 
wavelength dependence of extinction, they can be used to predict the extinction relationships for any photometric system by using 
synthetic photometry of artificially reddened energy distributions. An IDL procedure to produce the new curves at any desired value 
of R and over any wavelength range can be obtained from the author or via anonymous ftp at astrol.vill.edu. After logging in, 
change directories to pub/fitz/Extinction, and download the file "FMRCURVE.pro." Alternatively, this directory contains a series of 
compressed files named "FMRCURVEn.n.txt" (e.g., "FMRCURVE3.1.txt") which contain ASCII versions of the curves for various 
values of R ("n.n"). Two columns of data are contained in each file; the first contains wavelengths in A and the second contains the 
extinction curve in A(X)/E(B-V). There are 1099 sets of points in each file. 
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TABLE 1 

Uncertainties in Extinction Corrections 









Tm(X-V)/E(B-V) b 


Wavelength 


A" 1 


(for corrections using 


(for "ironing out" 


(A) 


(jxirT 1 ) 


average Galactic curve) 


the 2175 A bump) 


(1) 


(2) 


(3) 


(4) 


1000 


10.0 


2.5 


3.8 


1100 


9.1 


1.7 


2.8 


1200 


8.3 


1.3 


2.1 


1300 


7.7 


1.0 


1.7 


1400 


7.1 


0.86 


1.5 


1500 


6.7 


0.74 


1.3 


1600 


6.2 


0.66 


1.1 


1800 


5.6 


0.55 


1.0 


2000 


5.0 


0.55 


0.89 


2200 


4.5 


0.64 


0.78 


2400 


4.2 


0.43 


0.69 


2600 


3.8 


0.30 


0.62 


2800 


3.6 


0.23 


0.55 


3000 


3.3 


0.17 


0.48 


3500 


2.9 


0.06 


0.35 


4000 


2.5 


0.00 


0.25 



d k{\-V) = E(\-V)/E(B-V). These values of a k{X -v) should be inserted 
into Eq. 1 to compute the uncertainty in a corrected relative energy distribution 
m(X — V) when there is no information available about the true shape of the 
extinction law along the sightline and the average Galactic curve is adopted by 
default. See §3. in the text. 

b These values (multiplied by the derived E(B - V)) give the uncertainty in 
the final relative energy distribution m(X-V) when an object is dereddened by 
"ironing out" the 2175 A bump using the average Galactic extinction curve. 
See §3.1 in the text. 



Table 2 

Optic al/IR Extinction Ratios for R=3.1 



Extinction 


Observed 




Model Curve 


Ratio 


Value 


Reference 11 


Value 


(1) 


(2) 


(3) 


(4) 


A(M)/E(B-V) 


08-0 12 


1,2 


0.12 


A(L)/E(B-V) 


0.09-0.20 


1,2,3,4 


0.19 


A(K)/E(B-V) 


0.33-0.38 


2,3,4 


0.36 


A(H)/E(B-V) 


0.52-0.55 


1,2 


0.53 


A(J)/E(B-V) 


0.85-0.91 


1,2,3 


0.86 


A(I)/E(B-V) 


1.50 


3 


1.57 


A(R)/E(B-V) 


2.32 


3 


2.32 


A(V)/E(B-V) 


3.10 




3.10 


E(U-B)/E(B-V) 


0.70 + 0.05 xE(B-V) 


5 


0.69 + 0.04 xE(B-V) 


E(b-y)/E(B-V) 


0.74 


6 


0.74 


E(ml)/E(b-y) 


-0.32 


6 


-0.32 


E(cl)/E(b-y) 


0.20 


6 


0.17 


E(u-b)/E(b-y) 


1.5 


6 


1.54 



References: 1 = Rieke & Lebofsky 1985; 2 = Whittet 1988; 3 = Schultz & Wiemer 
1975; 4 = Savage & Mathis 1979; 5 = FitzGerald 1970; 6 = Crawford 1975 



Table 3 

Cubic Spline Anchor Points for R 



3.1 Curve 



Wavelength 


A" 1 


A(A) 


(A) 


(Mm" 1 ) 


E(B-V) 


(1) 


(2) 


(3) 


oo 


0.000 


0.000 


26500 


0.377 


0.265 


12200 


0.820 


0.829 


6000 


1.667 


2.688 


5470 


1.828 


3.055 


4670 


2.141 


3.806 


4110 


2.433 


4.315 


2700 


3.704 


6.265 


2600 


3.846 


6.591 
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Table 4 

R-Dependent Values of Optical Spline Anchor Points 



Wavelength 


A" 1 


A(A) 


(A) 


(Mm" 1 ) 


E(B-V) 


(1) 


(2) 


(3) 


6000 


1.667 


-0.426+1.0044x7? 


5470 


1.828 


-0.050+1.0016x7? 


4670 


2.141 


0.701 + 1.0016x7? 


4110 


2.433 


1 .208 + 1 .0032 x R - 0.00033 x 7? 2 



Table 5 

Comparison of Model IR Color Excesses With Observations 





E(B-V) 


E(V-J) 


E(V-H) 


E(V-K) 


E(V-L) 




Sightline 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


R 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


HD 210121 — observed 


0.38 


0.64 


0.71 


0.73 


0.78 




HD 210121 — model 




0.61 


0.70 


0.75 


0.79 


2.22 + 0.14 


HD 21483 — observed 


0.55 


1.07 


1.22 


1.32 


1.40 




HD 21483 — model 




1.08 


1.23 


1.31 


1.39 


2.69 + 0.11 


HD 167771 — observed 


0.44 


0.96 


1.14 


1.23 


1.25 




HD 167771 — model 




0.98 


1.13 


1.20 


1.27 


3.10 + 0.16 


HD 147889 — observed 


1.09 


2.98 


3.52 


3.90 


4.11 




HD 147889 — model 




3.11 


3.58 


3.81 


4.04 


3.92 + 0.17 


HD 37061 —observed 


0.54 


1.74 


2.05 


2.21 


2.34 




HD 37061 — model 




1.78 


2.06 


2.19 


2.32 


4.57 + 0.18 


HD 38087 — observed 


0.33 


1.10 


1.45 


1.58 


1.53 




HD 38087 — model 




1.21 


1.40 


1.49 


1.57 


5.08 + 0.58 


HD 36982 — observed 


0.34 


1.38 


1.61 


1.85 






HD 36982 — model 




1.43 


1.66 


1.76 




5.83 + 0.57 



NOTE. — Observed color excesses are derived using the intrinsic colors from Wegner 1994. Model color 
excesses are computed via synthetic photometry of artificially reddened stellar energy distributions using the 
IR extinction law derived in the Appendix and assuming the E(B — V) values from column 2 and the R values 
from column 7. The stated values of R yield the best fits to the observed IR color excesses and the uncertainties 
are based on the assumptions cfe(b-v) = 0.02 and UEiy-X) = 0.05. Kuracz ATLAS9 models with appropriate 
values of T e ff were used to represent the intrinsic stellar energy distributions. 
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FIG. 1 . — Normalized interstellar extinction curves from the far-IR through the UV. Several general features of the curves are noted. The solid and dotted curves 
are estimates for the case R = A(V)/E(B-V) = 3.1 derived in the Appendix of this paper and by Cardelli, Clayton, & Mathis 1989, respectively. The dashed curve 
shows the average Galactic UV extinction curve from Seaton 1979. 




3 4 5 6 7 8 9 

1/A (/im _1 ) 

FIG. 2. — Examples of 80 Galactic UV extinction curves derived from WE satellite observations. Analytical fits to the curves are shown, based on the work of 
Fitzpatrick & Massa 1990. The curves are taken from the Fitzpatrick & Massa catalog, with the addition of the lines-of-sight toward HD 210121 from Welty & 
Fowler 1992 and HD 62542 from Cardelli & Savage 1988. This figure demonstrates the enormous range of properties exhibited by UV extinction in the Milky Way. 
The dotted line, labeled "<r," shows the standard deviation of the sample scaled to the value <r(1500) = 0.74, as derived from ANS satellite data (see §3.1). 
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FIG. 3. — Analytical fitting functions for UV extinction curves from Fitzpatrick & Massa 1990. A normalized UV extinction curve (thick solid curve) can be 
represented by a combination of three functions: (1) a linear background component (thin dashed line), (2) a UV bump component (thin solid curve), and (3) a 
far-UV curvature component (thin dotted line). The linear background is parameterized by two tightly correlated coefficients (slope and intercept), the bump by 
three coefficients (strength, width, and central position), and the far-UV curvature by a single scale factor. Given the near invariance in bump central position and 
the correlation between the linear coefficients, most Galactic UV extinction curves can be reproduced to within the observational uncertainties with only four free 
parameters. 




FIG. 4. — Far-IR through UV extinction curves from Cardelli, Clayton, and Mathis 1989 (CCM). CCM found that extinction curves can be expressed approximately 
as a 1-parameter family that varies linearly with where R =A(V)/E(B-V) and has a mean value in the diffuse interstellar medium of 3.1. Examples of CCM's 
results are shown for four representative values of R, listed on the righthand side of the figure next to the corresponding curve. 
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FIG. 5. — Left Panel: Slope of the UV linear extinction component c2 plotted against R~ l (filled circles) for 31 sightlines from the Fitzpatrick & Massa 
1990 (FM) sample plus HD 210121 (at RT l = 0.45). The adopted linear relationship between these quantities is indicated with the dashed line and given by 
c2 = -0.824 +4.717/? -1 . Right Panel: Intercept of the UV linear extinction component cl plotted against the linear slope c2 (filled circles) for the full set of ~80 
extinction curves from the FM catalog. The adopted linear relationship between these parameters is indicated by the dashed line and given by cl = 2.030-3.007c2. 
The relationship between cl and c2 implicit in the CCM formula is shown by the dotted line. 
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FIG. 6. — A new estimate of the wavelength dependence of extinction in the IR/optical region for the case R = 3.1 (thick solid curve). Plus signs represent the 
extinction data from Bastiaansen 1992, normalized to R = 3.1; the dotted line shows the CCM curve for R = 3.1. The arbitrarily scaled profiles of the Johnson 
UBVRIJHKLM and Stromgren uvby filters are shown for comparison. The new curve was constrained to reproduce the broad- and intermediate-band filter-based 
extinction measurements listed in Table 2 and to fit the Bastiaansen data. For A > 2700 A (1/A < 3.7/im~' ) the curve is constructed as a cubic spline interpolation 
between the points marked by the filled symbols (see Table 3). At wavelengths shortward of 2700 A, the curve is computed using the FM fitting function with the 
coefficients given in the Appendix. 
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FIG. 7. — Examples of the R-dependent far-IR through UV extinction curves derived in this paper (thick solid and dashed curves). The corresponding values of R 
are listed on the righthand side of the figure beside the curves. For comparison, the results of CCM for the same four values of R are shown. 
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FIG. 8. — Left Panel: Observed-minus-predicted values of E(U -B)/E(B-V) plotted against the slope of the UV linear extinction component c2 (filled circles) 
for sightlnes from the FM sample. Model values of E(U - B) / E(B -V) vary over the range ~0.6 to ~0.8 for the observed range in c2. Right Panel: Observed 
minus predicted values of the Stromgren extinction ratio E(c\)/E(b-y) plotted against c2 (filled circles) for the FM sightlines with Stromgren data. Model values 
of E(c\)/E(f)-y) vary over the range ~0.0 to ~0.3 for the observed range in c2. 



